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"Warp drive" spacetimes are useful as "gedanken-experiments" that force us to confront the 
foundations of general relativity, and among other things, to precisely formulate the notion of 
"superluminal" communication. After carefully formulating the Alcubierre and Natario warp drive 
spacetimes, and verifying their non-perturbative violation of the classical energy conditions, we 
consider a more modest question and apply linearized gravity to the weak-field warp drive, testing 
the energy conditions to first and second order of the warp- bubble velocity, v. Since we take the warp- 
bubble velocity to be non-relativistic, u < c, we are not primarily interested in the "superluminal" 
features of the warp drive. Instead we focus on a secondary feature of the warp drive that has not 
previously been remarked upon — the warp drive (if it could be built) would be an example of a 
"reaction- less drive" . For both the Alcubierre and Natario warp drives we find that the occurrence of 
significant energy condition violations is not just a high-speed effect, but that the violations persist 
even at arbitrarily low speeds. 

A particularly interesting feature of this construction is that it is now meaningful to think of 
placing a finite mass spaceship at the centre of the warp bubble, and then see how the energy in 
the warp field compares with the mass-energy of the spaceship. There is no hope of doing this in 
Alcubierre's original version of the warp-field, since by definition the point in the centre of the warp 

C/-\ ■ bubble moves on a geodesic and is "massless" . That is, in Alcubierre's original formalism and in the 

Natario formalism the spaceship is always treated as a test particle, while in the linearized theory we 
can treat the spaceship as a finite mass object. For both the Alcubierre and Natario warp drives we 

\Q • find that even at low speeds the net (negative) energy stored in the warp fields must be a significant 

f*"^ ' fraction of the mass of the spaceship. 
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I. INTRODUCTION 



Alcubierre demonstrated, that it is theoretically possible, within the framework of general relativity, to attain 
arbitrarily large velocities [l|. In fact, numerous solutions to the Einstein field equations are now known that allow 
S^ "effective" superluminal travel J3, LJ HI Li B lj LlU LUl • Despite the use of the term superluminal, it is not "really" 

possible to travel faster than light, in any local sense. Providing a general global definition of superluminal travel is no 
trivial matter |12l Il3 | , but it is clear that the spacetimes that allow "effective" superluminal travel generically suffer 
from the severe drawback that they also involve significant negative energy densities. More precisely, superluminal 
effects are associated with the presence of exotic matter, that is, matter that violates the null energy condition 
[NEC]. In fact, superluminal spacetimes violate all the known energy conditions, and Ken Olum demonstrated that 
negative energy densities and superluminal travel are intimately related [l4| . Although most classical forms of matter 
are thought to obey the energy conditions, they are certainly violated by certain quantum fields |l5|. Additionally, 
certain classical systems (such as non-minimally coupled scalar fields) have been found that violate the null and the 
weak energy conditions |l6l H7I | . Finally, we mention that recent observations in cosmology strongly suggest that 
the cosmological fluid violates the strong energy condition [SEC] , and provides tantalizing hints that the NEC might 
possibly be violated in a classical regime [ji| UM, |2(j ■ 

For warp drive spacetimes, by using the "quantum inequality" deduced by Ford and Roman |2l| . it was soon 
verified that enormous amounts of energy are needed to sustain superluminal warp drive spacetimes |22l 1'2.']| . To 
reduce the enormous amounts of exotic matter needed in the superluminal warp drive, van den Broeck proposed a 
slight modification of the Alcubierre metric which considerably ameliorates the conditions of the solution 24] . It is 
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also interesting to note that, by using the "quantum inequality", enormous quantities of negative energy densities 
are needed to support the superluminal Krasnikov tube [S|. However, Gravel and Plante |25l l26| in a way similar in 
spirit to the van den Broeck analysis, showed that it is theoretically possible to lower significantly the mass of the 
Krasnikov tube. 

In counterpoint, in this article we shall be interested in applying linearized gravity to warp drive spacetimes, testing 
the energy conditions at first and second order of the warp-bubble velocity. We will take the bubble velocity to be 
non-relativistic, v < 1. Thus we are not focussing attention on the "superluminal" aspects of the warp bubble, such 
as the appearance of horizons [271 l28l |29j , and of closed timelike curves | 3('| , but rather on a secondary unremarked 
effect: The warp drive (if it can be realised in nature) appears to be an example of a "reaction-less drive" wherein the 
warp bubble moves by interacting with the geometry of spacetime instead of expending reaction mass. 

A particularly interesting aspect of this construction is that one may place a finite mass spaceship at the origin 
and consequently analyze how the warp field compares with the mass-energy of the spaceship. This is not possible 
in the usual finite-strength warp field, since by definition the point in the center of the warp bubble moves along a 
geodesic and is "massless" . That is, in the usual formalism the spaceship is always treated as a test particle, while in 
the linearized theory we can treat the spaceship as a finite mass object. 

Because we do not make any a priori assumptions as to the ultimate source of the energy condition violations, we 
will not use (or need) the quantum inequalities. This means that the restrictions we derive on warp drive spacetimes 
are more generic than those derived using the quantum inequalities — the restrictions derived in this article will hold 
regardless of whether the warp drive is assumed to be classical or quantum in its operation. 

We do not mean to suggest that such a "reaction-less drive" is achievable with current technology — indeed the 
analysis below will, even in the weak-field limit, place very stringent conditions on the warp bubble. These conditions 
are so stringent that it appears unlikely that the "warp drive" will ever prove technologically useful. The Alcubierre 
"warp drive" , and the related Natario "warp drive" , are likely to retain their status as useful "gedanken-experiments" 
— they are useful primarily as a theoretician's probe of the foundations of general relativity, and we wish to sound a 
strong cautionary note against over-enthusiastic mis-interpretation of the technological situation. 

II. WARP DRIVE BASICS 

Within the framework of general relativity, Alcubierre demonstrated that it is in principle possible to warp spacetime 
in a small bubble-like region, in such a way that the bubble may attain arbitrarily large velocities. Inspired by the 
inflationary phase of the early universe, the enormous speed of separation arises from the expansion of spacetime 
itself. The simplest model for hyper-fast travel is to create a local distortion of spacetime, producing an expansion 
behind the bubble, and an opposite contraction ahead of it. Natario's version of the warp drive dispensed with the 
need for expansion at the cost of introducing a slightly more complicated metric. 

The warp drive spacetime metric, in cartesian coordinates, is given by (with G = c = 1) 

ds 2 = —dt 2 + [dx — /3(x, y,z — zo(t)) dt] • [dx — fl(x, y, z — zo(t)) dt] . (1) 

In terms of the well-known ADM formalism this corresponds to a spacetime wherein space is flat, while the "lapse 
function" is identically unity, and the only non-trivial structure lies in the "shift vector" f3(t,x). Thus warp drive 
spacetimes can also be viewed as specific examples of "shift-only" spacetimes. The Alcubierre warp drive corresponds 
to taking the shift vector to always lie in the direction of motion 

(3(x, y,z- z (t)) = v(t) z f(x, y,z- z a (t)), (2) 

in which v(t) — dzo(t)/alt is the velocity of the warp bubble, moving along the positive z-axis, whereas in the Natario 
warp drive the shift vector is constrained by being divergence-free 

V-0(x,y,z) = O. (3) 

A. Alcubierre warp drive 

In the Alcubierre warp drive the spacetime metric is 

ds 2 = -dt 2 + dx 2 + dy 2 + [dz - v{t) f(x, y,z~ z {t)) dt} 2 . (4) 



The form function f(x, y, z) possesses the general features of having the value / = in the exterior and / = 1 in the 
interior of the bubble. The general class of form functions, f(x, y, z), chosen by Alcubierre was spherically symmetric: 
f(r) with r = \/x 2 + y 2 + z 2 . Then 

f(x,y,z-z (t)) = f(r(t)) with r(t) = {[(z - z {t)f + x 2 + y 2 } 1 ' 2 . (5) 

Whenever a more specific example is required we adopt 

tanh [a(r + R)] — tanh [a(r — R)] ,„. 

/(0 = 2tanh(ai?) ' (6) 

in which R > and a > are two arbitrary parameters. R is the "radius" of the warp-bubble, and a can be interpreted 
as being inversely proportional to the bubble wall thickness. If a is large, the form function rapidly approaches a top 
hat function, i.e., 

It can be shown that observers with the four velocity 

17" = (1,0,0,1;/), [/„ = (-!, 0,0,0). (8) 

move along geodesies, as their 4-acceleration is zero, i.e., a M = U v U^-^ = 0. These observers are called Eulerian 
observers. Their four- velocity is normal to the spatial hypersurface with t = const, and they are close analogues of 
the usual Eulerian observers of fluid mechanics, which move with the flow of the fluid. The spaceship, which in the 
original formulation is treated as a test particle which moves along the curve z = Zo(t), can easily be seen to always 
move along a timelike curve, regardless of the value of v(t). One can also verify that the proper time along this curve 
equals the coordinate time, by simply substituting z = Zo(t) in equation JIJ. This reduces to dr = di, taking into 
account dx = dy = and /(0) = 1. 

If we attempt to treat the spaceship as more than a test particle, we must confront the fact that by construction 
we have forced / = outside the warp bubble. [Consider, for instance, the explicit form function of equation JBJ in 
the limit r — > oo.] This implies that the spacetime geometry is asymptotically Minkowski space, and in particular the 
ADM mass (defined by taking the limit as one moves to spacelike infinity i°) is zero. That is, the ADM mass of the 
spaceship and the warp field generators must be exactly compensated by the ADM mass due to the stress-energy of 
the warp-ficld itself. Viewed in this light it is now patently obvious that there must be significant violations of the 
classical energy conditions (at least in the original version of the warp-drive spacetime) , and the interesting question 
becomes "Where are these energy condition violations localized?" . 

One of our tasks in the current article will be to see if we can first avoid this exact cancellation of the ADM 
mass, and second, to see if we can make qualitative and quantitative statements concerning the localization and 
"total amount" of energy condition violations. (A similar attempt at quan tification of the "total amount" of energy 
condition violation in traversable wormholes was recently presented in |3lll32j|.) 

Consider a spaceship placed within the Alcubierre warp bubble. The expansion of the volume elements, 8 = U^-^, 
is given by 9 = v (df/dz). Taking into account equation (JSJ, we have (for Alcubicrre's version of the warp bubble) 

e = v z-zoW (9) 

r dr 

The center of the perturbation corresponds to the spaceship's position zo(t). The volume elements are expanding 
behind the spaceship, and contracting in front of it. Appendix lAl contains a full calculation of all the orthonormal 
components of the Einstein tensor for the Alcubierre warp bubble. By using the Einstein field equation, G^ v — 8ir T^, 
we can make rather general statements regarding the nature of the stress energy required to support a warp bubble. 
The WEC states T^ v V^ V v > 0, in which V 11 is any timelike vector and T p „ is the stress-energy tensor. Its 
physical interpretation is that the local energy density is positive. By continuity it implies the NEC. In particular 
the WEC implies that T^ v U^ U v > where U^ is the four-velocity of the Eulerian observers discussed above. The 
calculations will be simplified using an orthonormal reference frame. Thus, from the results tabulated in Appendix lAl 
we verify that for the warp drive metric, the WEC is violated, i.e., 



- - v 2 
Tim U^U" = 



dx J \dy 



<0, (10) 



where T^p and U^ are, respectively, the stress energy tensor and timelike Eulerian four-velocity given in the orthonor- 
mal basis. Taking into account the Alcubierre form function ©, we have 



Tau U»U U = - 



1 v 2 (x 2 +y 2 ) 
32^ V 2 



dr 



-I 2 



<0. 



(11) 



By considering the Einstein tensor component, G t - t -, in an orthonormal basis (details given in Appendix |A"|) . and taking 
into account the Einstein field equation, we verify that the energy density of the warp drive spacetime is given by 



T tt = T »V jj ^ = rptt = j, = j, jjp, jjv ^ 



(12) 



and so equals the energy density measured by the Eulerian observers, that is, equation (|1 1|) . It is easy to verify that the 
energy density is distributed in a toroidal region around the z-axis, in the direction of travel of the warp bubble 23]. 
It is perhaps instructive to point out that the energy density for this class of spacetimes is nowhere positive. That 
the total ADM mass can nevertheless be zero is due to the intrinsic nonlinearity of the Einstein equations. 

We can (in analogy with the definitions in |3lLl32| ) quantify the "total amount" of energy condition violating matter 
in the warp bubble by defining 



M, 



warp 



Pwarp tl X 



T^ U^ U u d A x 



v 
32^ 



dr 



r 2 drd 2 fl = 

12 



dr 



r 2 dr. (13) 



This is emphatically not the total mass of the spacetime, but it characterizes how much (negative) energy one needs 
to localize in the walls of the warp bubble. For the specific shape function © we can estimate 



M, 



warp 



-v 2 R 2 a. 



(14) 



(The integral can be done exactly, but the exact result in terms of polylog functions is unhelpful.) Note that the 
energy requirements for the warp bubble scale quadratically with bubble velocity, quadratically with bubble size, and 
inversely as the thickness of the bubble wall. 

The NEC states that T^ W" k v > 0, where k^ is any arbitrary null vector and T M „ is the stress-energy tensor. Taking 
into account the Einstein tensor components presented in Appendix lAl and the Einstein field equation, G M „ = 87r T^, 
the NEC for a null vector oriented along the ±z directions takes the following form 



T h^ h u — 



dx 



dy 



In particular if we average over the ±z directions we have 

7: \Tfj, v k +g k +i + T^v k_~ k_~j 



8n \dx 2 dy 2 



8tt \dx) + \dy 



(15) 



(16) 



which is manifestly negative, and so the NEC is violated for all v. Furthermore, note that even if we do not average, 
the coefficient of the term linear in v must be nonzero somewhere in the spacetime. Then at low velocities this term 
will dominate and at low velocities the un-averaged NEC will be violated in either the +z or —z directions. 

To be a little more specific about how and where the NEC is violated consider the Alcubierre form function. We 

have 



1 tiu ^±z «-±z 



1 v 2 (x 2 +y 2 ) (d£ 
dr 



Sit 



± 



8- 



+ y 2 + 2(z-z {t)) 2 



dr 



y 2 d 2 f 



(17) 



The first term is manifestly negative everywhere throughout the space. As / decreases monotonically from the center 
of the warp bubble, where it takes the value of / = 1, to the exterior of the bubble, with / « 0, we verify that df /dr 
is negative in this domain. The term d 2 //dr 2 is also negative in this region, as / attains its maximum in the interior 
of the bubble wall. Thus, the term in square brackets unavoidably assumes a negative value in this range, resulting 
in the violation of the NEC. 

For a null vector oriented perpendicular to the direction of motion (for definiteness take k = ±x) the NEC takes 
the following form 



J-IIV K<±£ K±x 



8- 



9/ 
dy 



dz 



(1-/) 



d 2 f 
dz 2 



T 



( d 2 f 



8tt \dxdz 



(18) 



Again, note that the coefficient of the term linear in v must be nonzero somewhere in the spacetime. Then at low 
velocities this term will dominate, and at low velocities the NEC will be violated in one or other of the transverse 
directions. Upon considering the specific form of the spherically symmetric Alcubierre form function, we have 



~8tt 



y 2 + 2(z~z (t)) 2 fdf\ 2 M ,Jx 2 +y 2 d/ , (z - z (t)f d 2 / 



2 r 2 \ dr / V T 3 dr r 2 dr 2 

V2, 



v x(z-z (t)) fd 2 f Id/ 



l^v ^±x f"±x — 

T 8^ r 2 [~d7 2 ~rdr~)' ^ 

Again, the message to take from this is that localized NEC violations are ubiquitous and persist to arbitrarily low 
warp bubble velocities. 

Using the "volume integral quantifier" (as defined in [3l],[32]), we may estimate the "total amount" of averaged 
null energy condition violating matter in this spacetime, given by 

J V e ±i k v ± , d 3 x « J T„ v e ±i kl & d 3 x « -v 2 R 2 a^ M walp . (20) 

The key things to note here are that the net volume integral of the 0(v) term is zero, and that the net volume average 
of the NEC violations is approximately the same as the net volume average of the WEC violations, which are 0(v 2 ). 

B. Natario warp drive 

The alternative version of the warp bubble, due to Jose Natario demonstrates that the contraction/expansion 
referred to above is not always a feature of the warp drive [2|. In his construction a compression in the radial 
direction is exactly balanced by an expansion in the perpendicular direction, so that the shift vector is divergence-free 
V • (3 = 0. To verify this consider spherical coordinates (r, ip, </>) in the Euclidean 3-space (see Q] for details). The 
diagonal components of the extrinsic curvature are given by 

K rr = — 2v f ' cos ip and K vv = K^, = vf' cos <p , (21) 

respectively, where the prime denotes a derivative with respect to r. Taking into account the definition of the expansion 
of the volume element given by 9 = K l i , from equations (|2T)l . we simply have 9 = K rr + K vv + K^ = 0. This 
analysis provides a certain insight into the geometry of the spacetime. For instance, consider the front of the warp 
bubble, with cose/? > 0, such that K rr < (see [2| for details). This indicates a compression in the radial direction, 
which is compensated by an expansion, K vv + K^ = —K rr , in the perpendicular direction. Thus, this warp drive 
spacetime can be thought of as a warp bubble that pushes space aside and thereby "slides" through space. 

We now define (3 in the interior of the warp bubble to be v, and assume the radius of the warp bubble to be R. 
Furthermore we take to fall to zero over a distance of order 1/cr. Appendix iBl contains a full calculation of all the 
orthonormal components of the Einstein tensor for the Natario warp bubble, useful for analysis of energy condition 
violations in this class of spacetimes. 

For the WEC we use the results of Appendix to write 

V 17" W = i-G K = ~-±-Tr(K 2 ). (22) 

07T 1D7T 

This is manifestly negative and so the WEC is violated everywhere throughout the spacetime. 

We can again quantify the "total amount" of energy condition violating matter in the warp bubble by defining 

M warp = J pwarp d 3 x - J T^ U» U v d 3 x - - JL j Tr(K 2 ) d 3 x = -^1- I %,) /% J} d 3 x. (23) 

Now (3 in the interior of the warp bubble is by definition v, while the size of the warp bubble is taken to be R. 
Furthermore gradients of (3 in the bubble walls are of order va, and the thickness of the bubble walls is of order l/R. 
So we can again estimate (as for the Alcubierre warp bubble) 

Mwarp » -V 2 R 2 O. (24) 

With a bit more work we can verify that the NEC is violated throughout the spacetime. Consider a null vector 
pointing in the direction n, so that k^ = (l,fi J ). Then 

T MV k»k» = ±- [G & + 2n J G k + nW G i} ] (25) 



If the NEC is to hold, then this must be positive in both the n and — n directions, so that we must have 



\ {T, v e +h k\ h + T, v k\ kt h ) = ± [G U + WW Gts] > 



But now averaging over the x, y, and z directions we see that the NEC requires 



But from Appendix iBl we have 



G t ~ t - + ^GH = -2TV(K 2 



(26) 



(27) 



(28) 



which is manifestly negative, and so the NEC is violated everywhere. Note that K is 0(v) and so we again see that 
the NEC violations persist to arbitrarily low warp bubble velocity. 



III. LINEARIZED GRAVITY APPLIED TO THE WARP DRIVE 



Our goal now is to try to build a more realistic model of a warp drive spacctime where the warp bubble is 
interacting with a finite mass spaceship. To do so we first consider the linearized theory |33l l34l l35| applied to 
warp drive spacetimes, for non-relativistic velocities, v <C 1. A brief summary of linearized gravity is presented in 
Appendix IU1 In linearized theory, the spacetime metric is given by ds 2 = (rj^ + h^) dx^dx 1 ', with /i M „ <C 1 and 
Vnv — diag(— 1, 1, 1, 1). Taking into account equation (£IJ, it is an exact statement that /i M „ has the following matrix 
elements 



(K») = 



v 2 p -vf 





-vf 



(29) 



Now the results deduced from applying linearized theory are only accurate to first order in v. This is equivalent to 
neglecting the hoo — v 2 f 2 term in the matrix (/i M „), retaining only the first order terms in v. Thus, we have the 
following approximation 



(V) 



-vf 





-vf 



(30) 



The trace of h^ is identically null, i.e., h = h^^ = 0. Therefore, the trace reverse of h^, v , defined in equation l|C6|l . 
is given by h^ v — h^ v , i.e., equation (|30|l itself. 

The following relation will be useful in determining the linearized Einstein tensor, 



'():) 



dx a dxP 



d 2 f 
dx a dxP 



(31) 



where we consider the bubble velocity parameter to be constant, v = constant. For simplicity we shall use the original 
form function given by Alcubierre, equation (J5J, so that equation l|31|l takes the form 



'IKS 



dx a dxP 



dr dr d 2 / 
dx« dx? ch 5 " 



d 2 r d/ 



dx a dx@ dr 



(32) 



A. The weak energy condition (WEC) 



The linearized theory applied to Alcubierre's warp drive, for non-relativistic velocities (v <§; 1), is an immediate 
application of the linearized Einstein tensor, equation (|C7|) . and there is no need to impose the Lorenz gauge 36]. The 



interest of this exercise lies in the application of the WEC. In linearized theory the 4- velocity can be approximated 
by W = (1, 0, 0, 0), therefore the WEC reduces to 



and from equation l|C7|l . we have 



The respective terms are given by 



Tnu U^ U u — T o — 7T- Goo 

57T 



Goo — — 7: (^-oo^i M — h^^ 11 — 2 /io M ,o p ) 



hoo,^ = 0, 

hfjtv, = — 2/l03,03 j 

^0/i,0 A 



hrni.C! 11 — /io3,03 • 



(33) 



(34) 



(35) 
(36) 
(37) 



In general, if second time derivatives appear, they can be neglected because d/dt is of the same order as v d/dz, 
so that d^d^ = V 2 + 0{y 2 V 2 )). In fact, it is possible to prove that the Einstein tensor components Goo an d Gth 
do not contain second time derivatives of any generic h^ |33j. That is, only the six equations Gij — 8irTij, are 
true dynamical equations. In contrast, the equations Go M = SirTo^ are called constraint equations because they are 
relations among the initial data for the other six equations; which prevent one from freely choosing the initial data. 

Substituting equations (|33|) - (|3*7|l into equation (|31Jl . we have 



Thus, equation (|33|l is given by 



G 00 = O(v 2 ). 



T^ 17" U" = Too = 0(v 2 ) , 



(38) 



(39) 



and the WEC is identically "saturated" . Although in this approximation the WEC is not violated, it is on the verge of 
being so (to first order in v). This is compatible with the "exact" non-perturbative calculation previously performed. 



B. Negative energy density in boosted inertial frames 



Despite the fact that the observers, with U^ = (1,0,0,0), measure zero energy density [more precisely 0(v 2 )}, it 
can be shown that observers which move with any other arbitrary velocity, j3, along the positive z axis measure a 
negative energy density [at 0(v)]. That is, Tqq < 0. The f3 occurring here is completely independent of the shift vector 
P(x, y, z — zo(t)), and is also completely independent of the warp bubble velocity v. Consider a Lorentz transformation, 
x» = A^ x", with A^A"„=^ and A^ defined as 



(A">) 



with 7 = (1 — /3 2 ) x l 2 . The energy density measured by these observers is given by Too = ^a ^"d -^"" That is 



7 



7 /3" 

1 
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yp 


7 _ 



(40) 



T, 



no 



7 2 T 00 + 2 7 2 /3To3 + 7 2 /3 2 T 33 
7 2 /3 (2T 03 + /3T 3: 

7 2 /3 



(41) 



8tt 



2G 



I3G. 



■■',■■', 



(42) 



taking into account Too = 0, due to equation 139|) . The respective Einstein tensor components, G03 and G33, are 
given by 



G03 = -r (^03,11 + ^03,22) +0(v 2 ). 
G33 = Q(v 2 ). 



(43) 
(44) 



Finally, substituting these into equation (|42|l . to a first-order approximation in terms of v, we have 

7 2 /3 



T oo = — S^~ (ho3,u + ft-03,22) + 0(v ) 



Taking into account equation (|31|l . we have 



00 8tt ^<9x 2 <V 



0(^ 2 



Applying equation (|32f> . we have the following relations 

^03,11 = " 

^03,22 = - 



r 2 dr 2 \r r 3 J dr 
r 2 dr 2 \r r 3 J dr 



and substituting in equation l|45f) we therefore have 






j 2 /3v 
8tt 



d 2 / 
dr 2 



+ y 2 + 2(z-z (i)) i 



dr 



+ 0(v 2 ) 



(45) 
(46) 

(47) 
(48) 

(49) 



A number of general features can be extracted from the terms in square brackets, without specifying an explicit form 
of /. In particular, / decreases monotonically from its value at r = 0, / = 1, to / sw at r > R, so that df /dr is 
negative in this domain. The form function attains its maximum in the interior of the bubble wall, so that d 2 //dr 2 
is also negative in this region. Therefore there is a range of r in the immediate interior neighbourhood of the bubble 
wall that necessarily provides negative energy density, as seen by the observers considered above. Again we find that 
WEC violations persist to arbitrarily low warp bubble velocities. 

C. The null energy condition (NEC) 
The NEC states that T^ v k^ k v > 0, where fc M is a null vector. Considering k^ = (1, 0, 0, ±1), we have 



T h^h 1 ' 



1 

8^ 



(Gqo i 2Gq3 + G33) 



= T^- (/103, 11 



/J03,22 
2f Ffl 



8tt \dx 2 dy 2 



)+0(v 2 ) 
+ 0(v 2 ). 



(50) 



Whatever the value of the bracketed term, as long as it is nonzero the NEC will be violated (for small enough v) in 
either the +z or — z directions. This is compatible with the "exact" non-perturbative calculation previously performed. 
Using the general form function defined by Alcubierre, equation (|50|l reduces to 



1 n V K K ± 



Sn 



y 



d 2 / 

dr 2 



+ y 2 + 2{z-z (t)f 



dr 



+ 0(v 2 ) 



(51) 



Equation l|51|l is proportional to the energy density, Tgg, of equation (|49|l . We verify that the term in square brackets 
has a region for which it is negative, thus also violating the NEC in the immediate interior vicinity of the bubble wall. 



D. Linearized gravity applied to the spaceship 

The weak gravitational field of a static source, in particular of a spaceship, is given by the following metric 

ds 2 = -dt 2 + dx 2 + dy 2 + dz 2 - 2$(x, y, z) (dt 2 + dx 2 + dy 2 + dz 2 ) . (52) 



Applying the linearized theory with ds 2 = (r]^, v + h^) dcc M Ax u and h^ v <C 1, the matrix elements of hp, v are given by 



(K") 



-2$ 








-2$ 









-2$ 











-2$ 



(53) 



The trace is given by h = h^^ 



-4$. The elements of the trace reverse, /i„„, are the following 



(V) = 



-4$ 









Applying the linearized Einstein tensor, equation l|C7() . to the component, Goo, we have Poisson's equation 

V 2 $ = A-np, 



(54) 



(55) 



in which p is the mass density of the spaceship. This case is comparable to the Newtonian limit, which is valid when 
the gravitational fields are too weak to produce velocities near the speed of light, that is, |$| <C 1 and \v\ <C 1. For such 
cases general relativity makes the same predictions as Newtonian gravity. For fields that change due to the movement 
of the source, for instance along the z-axis with a velocity of v, we verify that d/dt is of the same order as v d/dz. It 
is easy to verify that the energy fluxes or the momentum densities, T$i are proportional to first order in v, and the 
stresses, ly, to second order. Therefore the stress-energy tensor T^ obeys the inequalities |Tbo| ^ \Toi\ 3> \Tij\. 



E. Spaceship immersed in the warp bubble 



Consider now a spaceship in the interior of an Alcubierre warp bubble, which is moving along the positive z axis 
with a non-relativistic constant velocity. That is, v <C 1. The metric is given by 

i2 



ds 2 = -di 2 + Ax 2 + Ay 2 + [Az - v f(x, y,z- vt) At]' 

-2$(x, y,z- vt) [At 2 + Ax 2 + Ay 2 + {Az ~ v f(x, y,z- vt) At) 2 



(56) 



If $ = 0, the metric H5fi|l reduces to the warp drive spacetime of equation (@J). If v = 0, we have the metric representing 
the gravitational field of a static source, equation l|52|). 

Applying the transformation, z' = Z — vt, the metric (|56|l takes the form 



ds 2 = -dt 2 + Ax 2 + dy 2 + [Az 1 - v {f(x, y, z) - 1} di] 2 

-2$(z, y, z') [At 2 + Ax 2 + Ay 2 + {Az' - v {f(x, y, z') - 1} di) 2 



(57) 



Note that the metric now looks "static" . Considering an observer in the interior of the bubble, co-moving with the 
spaceship, the metric (|57J) . with / = 1, reduces to 



As 2 



-At 2 



Ax 2 + Ay 2 + Az' 2 - 2$(x, y, z) [At 2 + Ax 2 + Ay 2 + Az' 2 ] 



so that inside the bubble we simply have the weak gravitational field of a static source, that of the spaceship. 
Outside the warp bubble, where / = 0, the metric l|57|l reduces to 



As 2 



-At 2 + Ax 2 + Ay 2 + [Az' + v Atf - 2$(x, y, z') [At 2 + Ax 2 + Ay 2 + (Az' + v At) 2 ] 



(58) 



(59) 



First order approximation 



Applying the linearized theory, keeping terms linear in v and <& but neglecting all superior order terms so that 
the approximation of equation (|C4|) is valid, the matrix elements, /i^„, of the metric (|56(l are given by the following 
approximation 



(V) 



-2$ -vf 

-2$ 

-2$ 

-vf -2$ 



(60) 
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The trace of h^ v is given by, h = h^ = —4$, and the trace-reversed elements, /i M „, take the following form 



iK") 



-4$ -vf 



-vf 



(61) 



Therefore, from the WEC, T tlu U^U v — T 00 , the trace-reversed elements that one needs to determine Goo [see 
equation 1341) ] are the following 



h ^ 



— hoofiO + /l00, 11 + ^00, 22 + ^00,33 j 
^00, 00 — 2/l 3,03 > 
ho/ifi 11 = ^^00,00 + ^03,03 ■ 

Substituting equations (|62|) - (|64[) into equation (|34[1 . we have 

Goo =2V 2 $ + 0(i) 2 ,w$,f 2 ). 
Thus, taking into account Poisson's equation, we verify that the WEC is given by 

T v JJ* t U v = p + 0(v 2 ,v<P,$ 2 ), 



(62) 

(63) 
(64) 

(65) 

(66) 



where p is now the ordinary energy density of the spaceship which is manifestly positive. In linearized theory, the 
total ADM mass of the space-time simply reduces to the mass of the space-ship, i.e., 

M ADM = fT 0Q d 3 x = J p d 3 x + 0{v 2 ,v<5>,<S> 2 ) = M ship + 0{v 2 ,v$,<S> 2 ) . (67) 

The good news for the warp drive aficionados is that the dominant term is manifestly positive. 
The NEC, with W = (1,0,0,±1), takes the form 



T h^V 



1 

8^ 



(Gqo =t 2Gq3 4- G33) . 



The Einstein tensor components, G03 and G33, are given by 



G03 - -- [hm^ -ho Mt3 -h 3li>0 



G 



33 



1 /— fXV — 

"2 y 1 ^., ~ 2/l 3/i,3 ; 



with the respective trace-reversed terms 



"■03,(U = ~^03, 00 + ^03,11 + ^-03,22 + ^03,33 
^0/i,3 = ~^00,30 + ^03,33 



'3/i,0 
- A 

Z 3^,3 



= -h. 
= -h. 



30,00 



30,30 



and the term h li u, ,iV is given by equation 1)63(1 . Substituting these into equations (I69() ~ ((70[) . we have 



G 



- = 3(0 + 0) ' ■**■■''■ 



G 33 = 0(v 2 ,v<S>,<5> 2 ). 

to first order in v and $, and neglecting the crossed terms v $. 
The null energy condition is thus given by 

v fd 2 f d 2 f 



877 V dx 



Vr = p ± — [-± + -± )+o(v 2 ,v$,&) 



dy 2 



(68) 

(69) 

(70) 



(71) 
(72) 
(73) 
(74) 



(75) 
(76) 



(77) 
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From this, one can deduce the existence of localized NEC violations even in the presence of a finite mass spaceship, 
and can also make deductions about the net volume-averaged NEC violations. First, note that for reasons of structural 
integrity one wants the spaceship itself to lie well inside the warp bubble, and not overlap with the walls of the warp 
bubble. But this means that the region where p ^ does not overlap with the region where the 0{v) contribution 
due to the warp field is non-zero. So regardless of how massive the spaceship itself is, there will be regions in the wall 
of the warp bubble where localized violations of NEC certainly occur. If we now look at the volume integral of the 
NEC, the 0(v) contributions integrate to zero and we have 



/t m „ fc£- k v ±i d 3 x = fp d 3 x + 0(v 2 ,v<I>,® 2 ) = M ship + 0{v 2 ,v®, $ 2 ) . 



(78) 



Consider a similar analysis for a null vector oriented perpendicular to the direction of motion, for instance, along 
the x— axis, so that k^ = (1, ±1, 0, 0). The NEC then takes the form 



T^WW = ±- (Goo ± 2G i + G„) . 
Taking into account equation (|C7|) . the respective Einstein tensor components are given by 

2 axoz 
Gii = 0{v 2 ,v<$>,<5> 2 ), 

and Goo is given by equation i|tj5fl . Thus, equation l|79|l takes the following form 

and considering the Alcubierre form function, we have 

57r r A \ dr z r dr I 



(79) 

(80) 
(81) 

(82) 
(83) 



As for the situation when we considered null vectors aligned with the direction of motion, for these transverse null 
vectors we find localized NEC violations in the walls of the warp bubble. We also find that the volume integral of the 
0(v) term is zero and that 



[T^k^k» ±i d 3 x= [ pd 3 x + 0(v 2 ,v<S>,$ 2 )=M ship + 0(v 2 ,v<I>,® 2 ). 



(84) 



The net result of this 0(v) calculation is that irrespective of the mass of the spaceship there will always be localized 
NEC violations in the wall of the warp bubble, and these localized NEC violations persist to arbitrarily low warp 
velocity. However at 0(v) the volume integral of the NEC violations is zero, and so we must look at higher order in 
v if we wish to deduce anything from the consideration of volume integrals to probe "net" violations of the NEC. 

Second order approximation 

Consider the approximation in which we keep the exact v dependence but linearize in the gravitational field of the 
spaceship $. The components of the Einstein tensor, relevant to determining the WEC and the NEC are 



Gtt = — 



' • : 



G ix 



yj .r:r 



O.r 



dx 



21 

dy 



dy 



2V 2 $ + 0($ 2 ) 

f o($ 2 ), 



1 (df\ 2 1 fdf 2 



dx 



dy 



Of 



(85) 

(86) 
(87) 



d 2 f 



dz I dz 2 



+ 0($ 2 



(89) 
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The WEC is given by 



T^ U fl U = p- 



v 
32^ 



dx 



dy 



or by taking into account the Alcubierre form function, we have 



0($ 2 ), 



0($ 2 ). 



327T r 2 y dr 

Once again, using the "volume integral quantifier" , we find the following estimate 

I Tp.o C/ A U° d 3 x = M ship -v 2 R 2 a + J 0($ 2 ) d 3 x , 

which we can recast as 

Mtotai = M ship + M warp + J 0($ 2 ) d 3 x , 



(90) 



(91) 



(92) 



(93) 



where M to tai is the net total integrated energy density, to the order of approximation considered, i.e., keeping the 
exact v dependence and linearizing in the gravitational field $. 

Now suppose we demand that the volume integral of the WEC at least be positive, then 



v 2 R 2 a < M shi . 



(94) 



This equation is effectively the quite reasonable condition that the net total energy stored in the warp field be less than 
the total mass-energy of the spaceship itself, which places a powerful constraint on the velocity of the warp bubble. 
Re-writing this in terms of the size of the spaceship i? s hi P and the thickness of the warp bubble walls A = 1/ct, we 
have 



, 2 < 



M s hi p -Rship A 



-Rshi 



R 2 



(95) 



For any reasonable spaceship this gives extremely low bounds on the warp bubble velocity. 
In a similar manner, the NEC, with fc M = (1, 0, 0, ±1), is given by 



T^k^k = p± 



jLf^i,^r\_^_ 



8ir \ dx 2 dy 2 J 8ir 



dx 



dy 



0($ 2 ). 



(96) 



Considering the "volume integral quantifier" , we verify that, as before, the exact solution in terms of polylogarithmic 
functions is unhelpful, although we may estimate that 



/ T AP fc A k° d 3 x = M ship - v 2 R 2 a + J 0($ 2 ) d 3 x , 



(97) 



which is [to order 0($ 2 )] the same integral we encountered when dealing with the WEC. This volume integrated NEC 
is now positive if 



v 2 R 2 a< M shi; 



(98) 



Finally, considering a null vector oriented perpendicularly to the direction of motion (for definiteness take k = ±x) , 
the NEC takes the following form 



1 fiu k k, — p 



8tt 



1 fdf 



2 \8y 



^V-(W)^ 

dz) [ J 'dz 2 



d 2 f 

8tt V dxdz 



0($ 2 



Once again, evaluating the "volume integral quantifier" , we have 



d/ 



dr 



T, k» W d 3 x = M ship - - -L r 2 dr + - (1-f) 2r-L+r 2 -L dr 



d/ , 2 d 2 / 



dr 



dr 2 



0($ 2 )d 3 a;, 



(99) 



(100) 
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which, as before, may be estimated as 

J Tp, fc A k° d 3 x w M ship - v 2 R 2 (7+j 0($ 2 ) d 3 x . (101) 

If we do not want the total NEC violations in the warp field to exceed the mass of the spaceship itself we must again 
demand 

v 2 R 2 a < M ship , (102) 

which places a strong constraint on the velocity of the warp bubble. 

IV. SUMMARY AND DISCUSSION 

In this article we have seen how the warp drive spacetimes (in particular, the Alcubierre and Natario warp drives) can 
be used as gedanken-experiments to probe the foundations of general relativity. Though they are useful toy models for 
theoretical investigations, as potential technology they are greatly lacking. We have verified that the non-perturbative 
exact solutions of the warp drive spacetimes necessarily violate the classical energy conditions, and continue to do so 
for arbitrarily low warp bubble velocity — thus the energy condition violations in this class of spacetimes is generic 
to the form of the geometry under consideration and is not simply a side-effect of the "supcrluminal" properties. 

Furthermore, by taking into account the notion of the "volume integral quantifier" , we have also verified that the 
"total amount" of energy condition violating matter in the warp bubble is negative. Using linearized theory, we have 
built a more realistic model of the warp drive spacetimc where the warp bubble interacts with a finite mass spaceship. 
We have tested and quantified the energy conditions to first and second order of the warp bubble velocity. By doing 
so we have been able to safely ignore the causality problems associated with "superluminal" motion, and so have 
focussed attention on a previously unremarked feature of the "warp drive" spacetime. If it is possible to realise even 
a weak-field warp drive in nature, such a spacetime appears to be an example of a "reaction-less drive" . That is, 
the warp bubble moves by interacting with the geometry of spacetime instead of expending reaction mass, and the 
spaceship (which in linearized theory can be treated as a finite mass object placed within the warp bubble), is simply 
carried along with it. We have verified that in this case, the "total amount" of energy condition violating matter 
(the "net" negative energy of the warp field) must be an appreciable fraction of the positive mass of the spaceship 
carried along by the warp bubble. This places an extremely stringent condition on the warp drive spacetime, namely, 
that for all conceivably interesting situations the bubble velocity should be absurdly low, and it therefore appears 
unlikely that, by using this analysis, the warp drive will ever prove to be technologically useful. Finally, we point 
out that any attempt at building up a "strong-field" warp drive starting from an approximately Minkowski spacetime 
will inevitably have to pass through a weak-field regime. Since the weak-field warp drives are already so tightly 
constrained, the analysis of this article implies additional difficulties for developing a "strong field" warp drive. In 
particular we wish to sound a strong cautionary note against over-enthusiastic mis-interpretation of the technological 
situation. 
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APPENDIX A: THE EINSTEIN TENSOR FOR THE ALCUBIERRE WARP DRIVE 

1. Generic form function 

The Einstein tensor is given by Gy, v = R^ — ^ g^u R- The Einstein tensor components of the warp drive spacetime, 
in cartesian coordinates, in an orthonormal basis, G/xo, with a generic form function, f(x, y, z — zo(t)), are given by 



G ff = -\v 2 



'it 



df\ 2 (dp 2 



dx ) \dy 



(Al) 



1 d 2 f 



G ** = -r^ ^ 
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G ty 



G^ 



^- T xx 

a* . 

^ xy 

(~rxz 



Gyz 






G z 



i d 2 f 

— v 

2 dydz ' 

2 V \dx 2 + dy 2 



dx 



dy 



dy 

dj_ 

dz 

d£ 
dz 



dy 



dx 

dx 

d£ 
dy 

2 



dj_ 
dz 



d 2 .f 



+ (1 - /) 



d 2 f 
dz 2 



2^ ^dxdz 
~ 2^^^ dydz 



' ( 1L 

dx 



dz 



(W)«4 



dz- 



dx 



dy 



(A3) 
(A4) 

(A5) 
(A6) 
(A7) 
(A8) 
(A9) 

(A10) 



2. Alcubierre form function 



Taking into account Alcubierre's choice of the form function, i.e., f(r) with r — y ' x 2 + y 2 + [z — zo(t)] 2 , the 
Einstein tensor components of the warp drive spacetime, in cartesian coordinates, in an orthonormal basis, G^, take 
the following form 



v 2 (x 2 + y 2 ) fdf 

« 4 r 2 y dr 

= v x(z-z (t)) fd 2 f Id/ 
t* 2 r 2 \dr 2 rdr 

= v y(z - z (t) ) (d 2 ! Id/ 
ty r,,2 



G iz - g 



dr 2 r dr 

+ y 2 d 2 / , x 2 +y 2 + 2(z-z (t)) 2 df 



dr 2 



v 2 x 2 - y 2 - 4(z - z (t)) 2 (d£ 2 

dr 





4 




C ~ - - 
'-' xy ~ 


v 2 xy 

T7 2 " 


m 


G X z - 


v 2 x(z 


-*>(*)) 


2 


r 2 


C ~~ - 

^•yz ~ 


v 2 y(z 


- zo(t)) 


2 


r 2 



dr 



+ v 2 (l-f) 



(z-z (t)) 2 d 2 / | x 2 +y 2 df 



2(^V-(l-/)f^--^ 

dr J \ dr 2 r dr 



d/ 
dr 



(1-/) 



W_ld/ 



v 2 y 2 - x 2 - 4(z - z (t)) 2 (dp ~ 



G yy ~ 4 r ^ 

3 2 (x 2 +y 2 ) fdf 

^zz — -V t, | —r- 

dr 



\ dr 2 r dr 

(z-z Q (t)) 2 d 2 / , x 2 +y 2 df 



%) +A1-/) 



dr 2 



r 3 dr 



(All) 
(A12) 
(A13) 
(A14) 
(A15) 
(A16) 
(A17) 
(A18) 

(A19) 
(A20) 
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3. Linearized Einstein tensor for the Alcubierre warp drive 

Linearizing in the warp velocity is very easy. The only nonzero components are 



1 d 2 f 



i oy 

— v 

2 dydz ' 



Giy = "V^, ( A22 ) 



1 fd 2 f d 2 f\ 

All other components are 0(v 2 ). 

For the Alcubierre form function f(r) we find that at linear order in warp bubble velocity 

_ _v x(z - z (t)) fd 2 f _ 1 d/\ 
W* - 2 r 2 ^ dr 2 rdr ) ' ^ AM ) 

_ _v y(z z (t)) fd 2 f 1 d A 



G fe - g 



a; 2 + y 2 d 2 / .t 2 + y 2 + 2(z - z (t)) 2 d/ 
r 2 dr 2 r 3 dr 



(A26) 



All other components are 0(v 2 ). 



APPENDIX B: THE EINSTEIN TENSOR FOR THE NATARIO WARP DRIVE 

The Einstein tensor for the Natario warp drive is relatively easy to calculate once one realizes that it can be dealt 
with within the framework of "shift-only" spacetimes. The extrinsic curvature of the constant-i time slices is 

Kij = /% ;J) = \ [dip, + djdi] (bi) 

and the no expansion condition implies 

K = Tr(K) = (B2) 

so the trace of extrinsic curvature is zero. The intrinsic curvature of the constant-t time slices is by construction zero 
(since they are flat). It is useful to introduce a vorticity tensor 

fiy = -/%,,] = \ [di0j - djfii] (B3) 

and then a tedious but straightforward computation (which follows a variant of the Gauss-Codazzi decomposition) 
yields 

Rijkl = K lk K i - KuK jk , (B4) 

Rtijk = -diiljk + 0i (KkiKij - KjiK ik ) , (B5) 

Rtitj = -dtKij + (Kn + flK).. - (K 2 ) tj - 0k0k,ij + k 0i (KkiKij - K jk K u ) . (B6) 

Here we have defined 

(KSI + flK).. = K ik n kj + n ik K k3 , (B7) 

and similarly 

(K 2 ) =K ik K kj . (B8) 



The computation reported here for the Natario warp drive is identical in spirit (and differs in only notation and 
minor technical details) to the computation of the "acoustic" Riemann tensor for the "acoustic geometry" reported 
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in [371 138| . Specifically, when sound waves propagate in a constant-density constant speed-of-sound background fluid 
flow the acoustic metric is identical in form to Natario's metric. 
The Einstein tensor is now (in an orthonormal basis) 



Git 



Gu 



Gr 3 



1. 



-Ti(K 2 






-6 i3 TT(K 2 ) - (K n + D.K) 



where d/di denotes the usual advective derivative 



di at ax j 



(B9) 
(BIO) 
(BH) 

(B12) 



In view of the fact that for the Natario spacetime all the time dependence comes from motion of the warp bubble, so 
that all fields have a space-time dependence of the form (3(x, y, z — zo{t)), we can deduce 



d 
dt 



dz dx l 



v(t)z - (3 



• V 



so that 



(B13) 

(B14) 

If we now linearize, then the only surviving term in the Einstein tensor is G^ = ^A/3j as all other terms are 0(v 2 ). 
A particularly simple special case of the Natario spacetime is obtained if the shift vector j3 is taken to be irrotational 
as well as being divergenceless J3g] In this special case f2 — > and also A/3j — ► 0, so that (in an orthonormal basis) 



Gv 3 



v(t) z - ■ VK i:J - -S tJ Tr(K 2 ) - {KD. + VtK) 



Ga 
Gn 



Gvj — 



1 



Tr(K 2 ) 



v{t)z-f3 ■VK ij --6 ij Tr{K 2 ). 



(B15) 
(B16) 
(B17) 



When linearizing this special case we find that all components of the Einstein tensor are 0(v 2 ). 



APPENDIX C: LINEARIZED GRAVITY 

For a weak gravitational field, linearized around flat Minkowski spacetime, we can write the spacetime metric 

asylum 



QfAV VptV ~r ft'liv ) 



(CI) 



with hav « 1. It is convenient to raise and lower tensor indices with rf v and rju V , respectively, rather than with 
g 1 *" and g^ v . We shall adopt this notational convention, keeping in mind that the tensor g^ v still denotes the inverse 
metric, and not r]^ a rf" g a p. It should also be noted that in linear approximation we have g^ v = rf v — W, since the 
composition of g^ v and g^ u differ from the identity operator only by terms quadratic in h^ v . 

The linearized Einstein equation can be obtained in the straightforward manner as follows. Adopting the form 
of equation i|Cl|) for the metric components, the resulting connection coefficients, when linearized in the metric 
perturbation h^, yield 



r M Q /3 = -ff" (h au ,p + hf3„ ia — h a p,v) 



(C2) 



A similar linearization of the Ricci tensor, R^ u — R a Mai ,, given as a function of the Christoffel symbols, i.e., 
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results in the following approximation 

where /i = ft, a a = rf-^ hp a is the trace of h a p. Applying a further contraction to discover the Ricci scalar, R = 
<7 M1/ .R M „ ~ rj^R^, the Einstein tensor to linear order is given by 

Equation l|C5(l can be simplified by defining the trace reverse of h a /3, given by 

/la/3 = h a /3 - -r/a/3 h , (C6) 

with h = h a = —h. Therefore, in terms of h a p, the linearized Einstein tensor reads 

1 iv 



G a p 2 



h a pJ L + VotpKv,^ - h a ^ - %,/ + O (h 2 a ^j , (C7) 

and the linearized Einstein equation is found to be 

- h a( 3.^ - riaph^^ + h afl>/ 3^ + hp^ = 167rT Q/3 . (C8) 

Without loss of generality, one can impose the gauge condition, i.e., 

TT,u = o, (C9) 

which is similar to the tensor analogue of the Lorenz gauge A^^ = of electromagnetism [36j . The Einstein equation 
then simplifies to: 

-ha(3,^ = 16nT a0 . (CIO) 

The gauge condition, equation (|C9|) . the field equation, equation l|C10|) . and the definition of the metric 

9/j.v = Viiv + hfiu = rjfj, v + h,^ — yViivh (Cll) 

are the fundamental equations of the linearized theory of gravity written in the Lorenz gauge 36] . 
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